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Abstract Here, a new species of bent-toed gecko, 
Cyrtodactylus phukhaensis sp. nov., is described from Doi 
Phu Kha, Nan province, Thailand based on molecular 
and morphological evidence. A phylogeny based on 
NADH dehydrogenase subunit 2 (ND2) and its flanking 
tRNAs places the new species in the chauquangensis group 
as a sister taxon to Cyrtodactylus wayakonei. The new 
species can be differentiated from other members of the 
chauquangensis group by having a unique combination 
of 7 or 8 supralabials, 23-28 longitudinal rows of dorsal 
tubercles, 8-10 infralabials, 9 femoral pores, 7 precloacal 
pores, 40-47 ventral scales, and a lack of bands crossing 
the temporal area. In addition, results from a chromosome 
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study of C. phukhaensis sp. nov. showed that the new 
species has a diploid chromosome number of 40 with a 
fundamental number of 46. The formula of the karyotype 
was as follows: 2n (40) = 2m + 4sm + 34t. Our findings 
suggest that further studies of Cyrtodactylus biodiversity in 
northern Thailand are needed. 


Keywords bent-toed gecko, Cyrtodactylus wayakonei, 
Cyrtodactylus chauquangensis, karyotype, phylogeny, 
taxonomy 


1. Introduction 


Cyrtodactylus, commonly known as bent-toed geckos, is the 
most diverse genus of geckonids with more than 320 described 
species (Uetz et al., 2021), many of which have been discovered 
ata high rate in the past decade. Cyrtodactylus inhabit the 
karstic region, a complex ecosystem known as a hotspot for 
vertebrate endemism and persistence (Grismer et al., 2021a). 


Therefore, many species of bent-toed geckos are endemic, 
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although they are widely distributed in India, Sri Lanka, the 
southern Himalaya region, Southeast Asia, Papua New Guinea, 
and some areas in northern Australia (Grismer et al., 2021a). 
Thus, changes to these unique and restricted habitats, such as 
those caused by wildfires, deforestation, and forest exploitation, 
can pose a high risk of extinction for the geckos. 

Presently, Cyrtodactylus is the most diverse genus in Gekkota 
with a wide distribution, but each species typically has a 
restricted or unique habitat (Grismer et al., 2020c). Studies of 
their biodiversity, however, are lacking because of their syntopic 
distribution with other geckonids. This distribution leads to 
cryptic or complex species, which make it difficult to determine 
whether unidentified specimens are new taxa or morphological 
variations. Several studies on Cyrtodactylus have confirmed the 
existence of cryptic species, which have since been reclassified 
as many new taxa using both morphological characterization 
and molecular tools (Grismer et al., 2012; Grismer et al., 2018; 
Nazarov et al., 2012; Sumontha et al., 2012). Therefore, the 
identification of Cyrtodactylus is still required, especially in 
some unexplored areas, to advance our understanding of the 
taxonomy, distribution, and biodiversity of Cyrtodact ylus spp. 
and facilitate their long-term conservation. 

According to the distribution of Cyrtodact ylus, Thailand is 
one area that harbors many Cyrtodactylus spp; these are mostly 
found in western and southern areas, which contain more 
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than half the species in Thailand (Uetz et al, 2021). However, 
a few Cyrtodactylus spp. have been found in the northern part 
of Thailand including Cyrtodactylus doisuthep, Cyrtodactylus 
dumnuii, Cyrtodactylus inthanon, Cyrtodactylus erythrops, 
Cyrtodactylus maelanoi, and Cyrtodactylus khelangensis, which 
were found in only three provinces (Chiang Mai, Mae Hong 
Son, and Lampang) (Bauer et al., 2009; Bauer et al, 2010; Grismer 
et al, 2020a; Kunya et al, 2014; Kunya et al, 2015; Pauwels et al., 
2014b). This finding is inconsistent with those for many lineages 
of Cyrtodactylus in the surrounding countries, e.g, Myanmar 
and Laos, which have continuous massif to those in northern 
Thailand with many limestone ecosystems that are suitable 
habitats. From several years of surveys in northern Thailand, 
an unknown Cyrtodactylus spp. was identified that differed 
from any previously described Cyrtodactylus spp. Therefore, 
this species of Cyrtodactylus collected from Nan province in 
northern Thailand was identified here as a new species based 
on its morphology, phylogeny, and karyotype. 


2. Materials and Methods 


2.1. Sample collection Samples were collected at Doi Phu 
Kha, Nan province, northern Thailand (Figure 1). Photographs 
of the specimens were taken to document their color pattern 
in life prior to euthanasia after which liver tissues were 


Figure 1 Distributions of the Cyrtodactylus chauquangensis group. Red star indicates the type locality of Cyrtodactylus phukhaensis sp. 
nov. in Nan Province, Thailand. Red dots represent the distributions of other Cyrtodactylus spp. including (1) C. zhenkangensis, (2) C. 
erythrops, (3) C. dumnuii, (4) C. doisuthep, (5) C. auribalteatus, (6) C. kunyai, (7) C. wayakonei, (8) C. vilaphongi, (9) C. spelaeus, (10) C. martini, (11) 
C. hekouensis, (12) C. sonlaensis, (13) C. otai, (14) C. puhuensis, (15) C. bobrovi, (16) C. huongsonensis, (17) C. cucphuongensis, (18) C. soni, and (19) 
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collected and stored in 9596 ethanol. The specimens were 
then preserved in 10% buffered formalin and stored in 70% 
ethanol. Type specimens were subsequently deposited at 
Agriculture University of Phayao, Phayao, Thailand (AUP). 
The geographical coordinates and elevations of the sampling 
locality were recorded using a Garmin eTrex 10 device (WGS84 
datum). 


2.2. Molecular study Total genomic DNA were extracted 
from liver tissues using a DNA extraction kit (RBC Bioscience, 
Taiwan, China). The NADH dehydrogenase subunit 2 (ND2) 
gene with partial flanking tRNA regions in mitochondrial 
DNA was amplified using PCR and the forward and reverse 
primers, L4437b (5-AAGCAGTTGGGCCCATACC-3) (Macey 
et al, 1997) and H5937b (5-AGGGTTCCGATATCTTTRTG-3) 
(Macey et al., 1999), respectively, in a total volume of 25 pL 


Table 1 List of specimens used in the phylogenetic analysis. 
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using 5x FIREPol® Master Mix (Solis BioDyne, Estonia). PCR 
reactions were completed under the following conditions: 
an initial denaturation at 95°C for 5 min and 35 cycles of 
denaturation at 95°C for 30 s, annealing at 50°C for 30 s, and 
extension at 72°C for 1 min followed by a final extension at 
72°C for 5 min. The PCR products were loaded onto 2% agarose 
gel, stained with GelRed® (Biotium, USA), and visualized under 
UV illuminator. The successfully amplified PCR products 
were then purified using a NucleoSpin® Gel and PCR Clean- 
up kit (Machery-Nagel, Germany) before being sequenced in 
both directions (A T G C, Bangkok, Thailand). The obtained 
sequences were checked manually and assembled using BioEdit 
version 7.0.5.3. The new sequences were deposited in the 
GenBank database under the accession numbers MN871845, 
MN871846, and MZ439914-MZ4399]7 (Table 1). 


Species Locality Voucher Accession No. Reference 

C. aequalis Kyaiktiyo Hill, Mon State, Myanmar LSUHC:12895 MF872275.1 Grismer et al. (2017) 

C. aequalis Kyaiktiyo Hill, Mon State, Myanmar LSUHC 14055 MN534918.1 — Grismer et al. (2020b) 

C. amphipetraeus Mae Sot District, Tak Province, Thailand AUP-00690 MT468912 Chomdej et al. (2021) 

C. auribalteatus Noen Maprang Eoo aac Province, AUP-01745 MZ439914 This study 

C. auribalteatus Ninen Megna Disma Prase mnll: Provine, AUP-01746 MZ439915 This study 

Thailand 

C. chauquangensis Quy Hop district, Nghe An province, Vietnam NA 2016.1 MT953475.1 Grismer et al. (2020c) 

C. dammathetensis Mawlamyine District, Mon State, Myanmar LSUHC:12862 MF872276.1 Grismer et al. (2017) 

C. dammathetensis Mawlamyine District, Mon State, Myanmar LSUHC:12864 MF872278.1 Grismer et al. (2017) 

C. doisuthep Mueang District, Chiang Mai Province, Thailand AUP-00777 MT497801 Chomde;j et al. (2021) 

C. dumnuii Chiang Dao District, Chiang Mai Province, Thailand AUP-00769 MT497802 Chomdej et al. (2021) 

C. dumnuii Chiang Dao District, Chiang Mai Province, Thailand AUP-00770 M'T497803 Chomdej et al. (2021) 

C. erythrops Pang Mapha.D i UNE Son Province; AUP-00771 MT497806 ^ Chomdejetal. (2021) 

C. huongsonensis My Duc District, Hanoi, Vietnam IEBR A.2011.3A MT953481.1 Grismer et al. (2020c) 

C. inthanon Chom Thong District, Chiang Mai Province, Thailand —— AUP-00156 MT497800 Chomdej et al. (2021) 

C. kunyai Nong Hin District, Loei Province, Thailand AUP-01747 MZ439916 This study 

C. kunyai Nong Hin District, Loei Province, Thailand AUP-01748 MZ439917 This study 

C. phukhaensis sp. nov. Doi Phu Kha, Pua District, Nan Province, Thailand AUP-01823 MNS871845 This study 

C. phukhaensis sp. nov. Doi Phu Kha, Pua District, Nan Province, Thailand AUP-01824 MNS$871846 This study 

C. puhuensis Thanh Hoa Province, Vietnam ND 01.15 MT953489.1 Grismer et al. (2020c) 

C. sinyineensis Hpa-an District, Kayin State, Myanmar LSUHC:12835 MF872354.1 Grismer et al. (2017) 

C. sinyineensis Hpa-an District, Kayin State, Myanmar LSUHC:12836 MF872355.1 Grismer et al. (2017) 

C. sinyineensis Hpa-an District, Kayin State, Myanmar LSUHC:12837 MF872356.1 Grismer et al. (2017) 

C. soni Gia Vien District, Ninh Binh Province, Vietnam TEBR R.2016.4 MT953491.1 Grismer et al. (2020c) 

C. sonlaensis Phu Yen District, Son La Province, Vietnam IEBR A.2017.1 MT953492.1 Grismer et al. (2020c) 

C. vilaphongi Leang Prabang S d Prabangromcs TERRA 201G103 | MIT953497.11 | Gristuer eh al (Q020c) 

C. wayakonei Vieng Phoulkia ddp eds Nam Tha Province,  2EMK91016 MT953498.1  Grismer et al. (2020c) 

Outgroup 

Dixonius siamensis Thong Pha X hum oe Province,  AUP-01724 MT468896.1 ^ Chomdej et al. (2021) 

Gekko gecko Mimbu Township, Magway Division, Myanmar CAS 213628 JNO19053.1 Rósler et al. (2011) 
Thong Pha Phum District, Kanchanaburi Province, g 

Gekko kaengkrachanense AUP-01710 MT468895.1 Chomdej et al. (2021) 


Thailand 
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2.3. Phylogenetic analysis Phylogenetic relationships were 
analyzed using both maximum likelihood (ML) and Bayesian 
inference (BI) methods. In our analysis, we used the new 
sequences obtained in this study, additional sequences (taken 
from GenBank) of currently known Cyrtodactylus distributed 
close to the sampling site, and outgroups (Gekko gecko, Gekko 
kaengkrachanense, and Dixonius siamensis; Table 1). The sequences 
were aligned using ClustalW, and then manually checked and 
trimmed in MEGA X version 10.2.2 (Kumar et al., 2018). The 
best fit substitution models for ML and BI were estimated 
using jModelTest 2.1.7 (Darriba et al., 2012; Posada, 2008), which 
supported GTR+I+G for the ND2 gene and flanking tRNA 
regions as suggested by the Akaike information criterion. 
MrBayes version 3.2.7 (Ronquist et al., 2012) was used for BI 
analysis with a Markov chain Monte Carlo method involving 
10 000 000 generations sampled every 1 000 generations. The 
first 2 500 trees were discarded as burn in and the BI tree was 
visualized in ITOL version 6.11 (http;//itoLembl.de). ML analysis 
was performed using MEGA X version 102.2 (Kumar et al., 
2018) with 1 000 bootstrap replicates. The confidence of the tree 
topology was assessed using nodes that had Bayesian posterior 
probabilities (BPP) 20.95, while bootstrap support (BS) values 
>95 were considered as strong support. Nodes with BPP and 
BS values of 0990-094 were considered well supported. Mean 
uncorrected pairwise sequence divergence among and between 
species was calculated using MEGA X version 10.0.5 (Kumar et 
al, 2018). 


2.4. Morphological study Measurements were taken with 
digital calipers to the nearest 0.1 mm. Bilateral scale counts 
were given as left/right. We followed the methodology of 
measurements and meristic counts described by Grismer et 
al. (2017). The following terms were abbreviated as shown 
in parentheses: snout-vent length (SVL), tail length (TL), tail 
width (TW), forearm length (FL), tibia length (TBL), axilla-to- 
groin length (AG), head length (HL), head width (HW), head 
depth (HD), eye diameter (ED), eye to snout distance (ES), eye- 
to-nostril distance (EN), interorbital distance (IO), ear diameter 
(EL), internarial distance (IN), supralabials (SL), infralabials 
(IL), ventral scales (VS), subdigital lamellae on the fourth toe 
(ET4), femoral pores (FP), precloacal pores (PP), postprecloacal 
scales rows (PPS), body bands (BB), light bands on the tail (LB), 
and dark bands on the tail (DB). Morphological comparisons 
and analyses were based on specimen examination and data 
obtained from the existing literature (Bauer et al., 2009; Bauer et 
al, 2010; Kunya et al, 2014; Le et al, 2016; Liu and Rao, 2021; Luu 
et al, 2011; Nguyen et al., 2014; Nguyen et al, 2017; Nguyen et al, 
2010; Pauwels et al, 2014a; Quang et al, 2007; Schneider et al., 
2014; Sumontha et al, 2010; Tri, 2011; Zhang et al., 2021). 


2.5. Chromosome preparation and staining Metaphase 
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chromosomes were prepared using the squash technique 
adapted from Pradit et al. (2018). An intraperitoneal injection 
of colchicine solution (1 mg/mL) was administered at 15 ug/g 
body weight for 4 h. Subsequently, mitotic and meiotic cells 
were collected from the bone marrow of the fore (radius, ulna, 
and humerus bones) and hind limbs (fibula, tibia, and femur 
bones) as well as the testes. These cells were incubated in 0.075 M 
KCI solution for 30 min and then fixed in 3:1 methanol/acetic 
acid solution for three times. The cell spread was prepared 
by dropping 1 or 2 drops of cell suspension onto a precooled 
slide and allowing the suspension to air dry. The chromosome 
spreads were stained using 10% Giemsa solution for 15 min, and 
the well-spread chromosomes in metaphase were observed and 
photographed under a light microscope (Olympus Corporation, 
Tokyo, Japan). Karyotypes and ideograms were analyzed and 
constructed using KaryoType software (Altinordu et al, 2016). 


3. Results 


3.1. Molecular relationship The topology using the ND2 
fragment and flanking tRNA regions derived from Bayesian 
analysis was highly similar to that from ML analysis. The 
genealogical relationships of Cyrtodact ylus in this study were 
revealed to contain two clades (Figure 2); one consisted of 
Cyrtodact ylus spp. from northern and western Thailand, 
including nearby regions in Mon and Karen, Myanmar, 
which included members of the Cyrtodactylus sinyineensis group 
corresponding to a previous report by Grismer et al. (2021b); the 
other contained Cyrtodactylus from northern and northeastern 
Thailand, including northern Laos and Vietnam. The newly 
discovered species, C yrtodact ylus phukhaensis sp. nov. from 
Nan, northern Thailand, was nested in this clade with a clearly 
divergent lineage. The newly identified species represented a 
close sister taxon to a known species, C. wayakonei, with strong 
levels of nodal support (100/100 for BPP/BS values). This clade 
comprised C yrtodact ylus members of the C. chauquangensis 
group, as reported by Grismer et al. (2021b). 


32. Genetic divergence The uncorrected p-distances in the 
ND2 gene with flanking tRNA regions among and within the 
Cyrtodactylus spp. examined in this study are shown in Table 2. 
The observed interspecific distance within Cyrtodactylus varied 
from 5.2296 between C. phukhaensis sp. nov. and C. wayakonei 
to 33.71% between Cyrtodactylus sonlaensis and C. sin yineensis. 
The genetic difference among C. phukhaensis sp. nov. and other 
Cyrtodactylus species was 5.22%-24.73%. The species most closely 
related to C. phukhaensis sp. nov. were C. wayakonei (5.22%) and 
C. dumnuii (13.76%). Intraspecific genetic variations in the ND2 
gene with flanking tRNA regions was observed to be very 
low (0.00% in C. phukhaensis sp. nov., Cyrtodact ylus auribalteatus, 
Cyrtodactylus kun yai, Cyrtodactylus aequalis, C. sin yineensis, and 
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Figure 2 Bayesian inference phylogenetic tree of the mitochondrial ND2 gene with flanking tRNAs from the species of Cyrtodactylus in 
the chauquangensis and sinyineensis groups. Support values at the nodes represented the Bayesian posterior probability and maximum 


likelihood bootstrap support values from the same dataset. 


Cyrtodactylus dammathetensis to 0.01% in C. dumnuii). The newly 
identified species, C. phukhaensis sp. nov., was clearly distinct 


from other known species of Cyrtodactylus. 
3.3. Morphological analyses 


Taxonomic accounts 

Cyrtodactylus phukhaensis sp. nov. Chomdej, Pradit, 
Pawangkhanant, Naiduangchan, and Suwannapoom 

Phukha bent-toed gecko (Figures 3-7) 

Holotype AUP-01823 (Figures 3 and 4), adult male from Doi 
Phu Kha, Pua District, Nan Province, Thailand, (19910.40'N; 
101906.44'E; 1 700 m as.l.), collected on October 25, 2018 by 
Parinya PAWANGKHANANT and Mali NAIDUANGCHAN. 


Paratypes AUP-01824 and AUP-01825 (Figure 5), two adult 
females; collected between October 25 and 26, 2018 by Parinya 
PAWANGKHANANT and Mali NAIDUANGCHAN from 
Doi Phu Kha, Pua District, Nan Province, Thailand (19910.40'N; 
101°06.44E; 1 700 m asl). 

Etymology The specific epithet, phukhaensis, is derived from 
the place name of Doi Phu Kha where the type series was 
collected. 

Diagnosis Cyrtodactylus phukhaensis sp. nov. differs from 
all other species in the C. chauguangensis group by having a 
combination of 8 or 9 supralabials; 8-10 infralabials; 36-44 


paravertebral tubercles; 23-28 longitudinal rows of dorsal 
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Figure 3 Dorsal view of live Cyrtodactylus phukhaensis sp. nov. 
from Doi Phu Kha, Pua and Bo Kluea District, Nan Province, 


Figure 4 Holotype of adult male Cyrtodactylus phukhaensis sp. 


nov. (AUP-01823). (A) Lateral view of head; (B) dorsal view 
of top of head; (C) precloacal region; (D) precloacal-femoral 
series showing scalation and pore arrangements; (E) gular 
region showing mental, postmental, chin scale arrangement, 


and pair of chinshields under dentary; (F) plantar view of feet 


showing subdigital lamellae under left hindfeet toes. 


on the fourth toe, 11 or 12 unmodified subdigital lamellae on 
the fourth toe, and 18 or 19 total subdigital lamellae on the 


tubercles; 40-47 ventral scales; 7 expanded subdigital lamellae 
pore-bearing precloacal scales in males; 3 rows of enlarged 


fourth toe; 27-34 enlarged femoral scales; 10-12 pore-bearing 
femoral scales in males; 8-11 enlarged precloacal scales; 7 


postprecloacal scales; 4-7 broken to hour glass-shaped dorsal 
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Figure 5 Dorsal view of paratypes of Cyrtodactylus phukhaensis 
Sp. nov. in preservative. 


Figure 6 Habitat of Cyrtodactylus phukhaensis sp. nov. at the 
Doi Phu Kha, Nan Province, Thailand, i.e., the type locality. 


body bands; 4-8 light-colored caudal bands (n = 2); 3-8 dark- 
colored caudal bands (n = 2); raised and strongly keeled dorsal 
tubercles that extend beyond the base of the tail; enlarged 
femoral and precloacal scales that are nearly the same size and 
continuous; pore-bearing femoral and precloacal scales that 
are not continuous; medial subcaudals 2-3 times wider than 
long and extending onto the lateral side of the tail; green irises; 
a nuchal loop lacking an anterior azygous notch and bearing 
a jagged posterior border; dorsal bands bearing paravertebral 
elements that are generally equal in width to the interspaces, 
bear lightened centers, and are edged with white tubercles; dark 
markings in the dorsal interspaces; light caudal bands in adults 
bearing dark-colored markings; light-colored caudal bands that 
do not encircle the tail, and a mature regenerated tail that is 
not spotted (Table 3). 
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Figure 7 The specimen observed on a large boulder near the 
stream in the habitat. 


Description of holotype Adult male, SVL = 778 mm (Figure 
4); head moderate in length (HL/SVL = 028), wide (HW/HL = 
0.66) and flat (HD/HL = 0.28), distinct from neck, triangular in 
dorsal profile; lores inflated, prefrontal region concave, canthus 
rostralis rounded; snout, elongated (ES/HL - 0.39), rounded 
in dorsal profile, broad in lateral profile; eye, large (ED/HL 
= 0.25); ear opening, oval (EL/HL = 0.10); eye-to-ear distance, 
greater than eye diameter; rostral, rectangular, partially divided 
dorsally, bordered posteriorly by supranasals and one internasal, 
laterally by first supralabials; external nares, bordered 
anteriorly by rostral, dorsally by supranasals, posteriorly by 
smaller postnasals, and ventrally by first supralabials; 9 (R,L) 
rectangular supralabials, extending to below eye midpoint; 7 
(R,L) infralabials, tapering posteriorly to commissure of jaw; 
scales of rostrum and lores, slightly raised, larger than granular 
scales on top of head and occiput; scales on top of head and 
occiput, intermixed with small tubercles; dorsal superciliaries, 
weakly pointed and directed posteriorly; mental triangular, 
bordered laterally by first infralabials and posteriorly by large 
left and right trapezoidal postmentals, which contact medially 
for 50% of their length posterior to mental; one row of enlarged 
chinshields, outermost row bordering first five infralabials; 


gular and throat scales, granular, grading posteriorly into 
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Table 3 Measurements (mm), meristic data, and color pattern for the type series of Cyrtodactylus phukhaensis sp. nov. 


AUP-01823 AUP-01824 AUP-01825 

Sex M F E 
SVL 77.8 82.07 82.8 
iL; 80.99 75.81 53.74 
TW 5.42 5:35 6.75 
EI 2217/5] 22.38 24.71 
TBL 8.46 11.2 10.01 
AG 3152 S 35.81 
HL 21.45 23.73 23.31 
HW 13.64 14.94 16.69 
HD 6.13 6.38 6.9 
ED 4.32 4.02 5.58 
ES 9.16 10.11 10.21 
EN 6.67 T 7.98 
IO 5.99 6.22 6.5 
EL 1.96 15772) 11208) 
IN 2.32 2.64 1.86 
SIL, 8 9 8 
IL 8 10 9 
PV 39 44 36 
LT 28 23 24 
VS 47 40 40 
ET4 12 12 12 
BD 9 D 9 
PP 7 7 7 
PPS 3 3 3 
BB 5 5 5 
LB 8 4 - 
DB 8 3 - 
Tubercles on at least anterior 1/3 of tail no no no 
Transverse subcaudal 2-3 times wider than long yes yes no 
Dorsal pattern with paravertebral elements yes yes yes 
Paravertebral elements not in contact yes yes yes 


larger, subimbricate pectoral and ventral scales. Body relatively 
short (AG/SVL = 0.42) with well-defined ventrolateral folds; 
dorsal scales, small, raised, and interspersed with large, raised, 
semi-regularly arranged, strongly keeled tubercles; tubercles, 
extending from top of head to base of tail just beyond the 
postcloacal swelling; tubercles on nape, smaller than those on the 
body; 35 paravertebral tubercles; approximately 18 longitudinal 
rows of dorsal tubercles; 28 flat, subimbricate, ventral scales 
larger than dorsal scales; 11 enlarged precloacal scales; 8 pore- 
bearing precloacal scales, separated on the midline by 1 poreless 
scale; 3 rows of large, postprecloacal scales; and no deep 
precloacal groove or depression. Forelimbs moderate in stature, 
relatively short (FL/SVL = 0.17); forearm scales, slightly raised, 
juxtaposed, larger than those on body, intermixed with large 
tubercles; palmar scales, slightly raised, juxtaposed; digits, well- 
developed, relatively long, inflected at basal interphalangeal 
joints; digits narrow distal to inflections; widened proximal 
subdigital lamellae extending onto palm; slight webbing at base 


of digit; claws, well-developed, sheathed by a dorsal and ventral 
scale at base; hind limbs more robust than forelimbs, moderate 
in length (TBL/SVL = 0.15), covered dorsally by small, raised, 
juxtaposed scales intermixed with large pointed tubercles and 
bearing flat, slightly larger imbricate scales anteriorly; ventral 
femoral scales flat, imbricate, much larger than dorsals; one 
row of 13 or 14 (R,L) enlarged femoral scales and enlarged 
precloacal scales, continuous; enlarged femoral scales, nearly 
equal in size; small postfemoral scales form an abrupt union 
with larger, flat, ventral scales on posteroventral margin of 
thigh; 6 (R,L) pore-bearing femoral scales, not continuous with 
pore-bearing precloacal scales; subtibial scales, flat, imbricate; 
plantar scales, raised; digits, relatively long, well-developed, 
inflected at basal, interphalangeal joints; 7 (R,L) transversely 
expanded subdigital lamellae on fourth toe proximal to joint 
inflection, not extended to sole; 11 (R,L) unmodified subdigital 
lamellae, distal to inflection; claws well-developed, base of 


claw sheathed by a dorsal and ventral scale. Tail complete, 
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original, gracile in proportions, 115.1 mm in length, 9.9 mm 
in width at base, tapering to a point (TL/SVL - 142); dorsal 
scales of tail flat, forming indistinct whorls; median row of 
transversely expanded subcaudal scales, three times as wide as 
long, extending onto lateral subcaudal region; three enlarged 
postcloacal tubercles at base of tail on hemipenal swellings; 
postcloacal scales large and flat. 

Variation Paratypes are similar to holotypes in terms of the 
morphological characteristics, measurements, and scalations 
provided in Table 3. Despite a limited number of specimens, 
few variations were found in the three geckos. Labrials varied 
in relation to supralabials 8-9 and infralabials 8-10 (Figure 5). 
Distribution Cyrtodactylus phukhaensis sp. nov. is known 
only from the type locality at Doi Phu Kha, Pua and Bo Kluea 
District, Nan Province, northern Thailand. 

Ecology and natural history The type series was collected 
at night between 19:00 and 23:00 in the Montane forest at 
elevations of 1 250-1715 m (Figure 6). Specimens were found on 
large rocks and large trees about 10-2.5 m above the stream 
and uphill. In addition, one specimen was observed on a large 
boulder near a stream (Figure 7). 

Only one adult male (holotype) was collected at night from 

the large boulder near the stream, which was covered with 
lower herbaceous plants. The male was sitting upside down on 
the boulder 15 m above the stream. Another specimen observed 
at the summit of the national park was observed on a large 
Sapindaceae (Acer sp.) tree. One of the paratypes (AUP-01825) 
was observed to predate an unknown orthopteran. Other 
herpetofauna species were observed as follows: Pseudocalotes 
microlepis (Agamidae), Dopasia gracilis (Anguiidae), Elaphe 
taeniura cf. mocquardi (Colubridae), Ovophis monticola (Viperidae), 
and Rhacophorus kio (Rhacophoridae). 
Comparisons Species of Cyrtodactylus are known to be 
endemic to their narrow distributions because of geoisolation; 
thus, the comparisons in this study were mainly between the 
new species and those from the C. chauquangensis group in its 
neighborhood, typically Cyrtodactylus spp. from northern and 
northeastern Thailand, northern Laos, and Vietnam (Figure 1). 

Based on phylogeny, the new species of Cyrtodactylus is a 
sister to C. wayakonei (p-distance = 5.22%, Table 2 and Figure 
2). Morphologically, C. phukhaensis sp. nov. differed from C. 
wa yakonei in the following ways: the number of supralabials (C. 
phukhaensis sp. nov. 8 or 9 vs. C. wayakonei 7 or 8); the number 
of body tubercle rows (C. phukhaensis sp. nov. 23-28 vs. C 
wa yakonei 17-19); the condition of femoral pores (distinct on C. 
phukhaensis sp. nov. vs. absent on C. wa yakonei); and the number 
of precloacal pores (C. phukhaensis sp. nov. 7 vs. C. wa yakonei 
6—8). Moreover, C. phukhaensis sp. nov. can also be distinguished 
from other neighboring C. chauquangensis group species by 
having the following: 8 or 9 supralabials (vs. C. auribalteatus 
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10 or 11, C. chauquangensis 9 or 10, C. doisuthep 10-12, and C. 
huongsonensis 10-13); 8-10 infrala bials (vs. C. chauquangensis 9-11 
and C. huongsonensis 10-12); the absence of a band crossing the 
temporal area (present on C. auribalteatus, C. chauquangensis, C. 
doisuthep, C. erythrops, C. huongensis, C. kunyai, C. puhuensis, C. 
sonlaensis, and C. vilaphongi); 23-28 rows of trunk tubercles (vs. 
C. doisuthep 19-20, C. dumnuii 18-22, C. huongsonensis 14-16, C. soni 
10-13, C. sonlaensis 13-15, and C. vilaphongi 15 or 16); 5 narrow 
body bands (n = 3) (vs. 4 bands on C. auribalteatus, irregular 
broad bands on C. chauquangensis, 6 or 7 bands on C. doisuthep, 
and 6 or 7 bands on C. dumnuii); 40-47 ventral scales (vs. C. 
auribalteatus 38—40, C. chauquangensis 36-38, C. doisuthep 29-35, C. 
erythrops 28, C. kunyai 34, C. puhuensis 36, and C. vilaphongi 34— 
36); 9 femoral pores (vs. C. auribalteatus 4 or 5, C. chauquangensis 
O, C. doisuthep O, C. dumnuii 6, C. erythrops 9 or 10 on each side of 
thigh, C. huongsonensis 15-17, C. kun yai 10-12, C. puhuensis 0, C. 
soni 12-16, and C. sonlaensis 28-30); and 7 precloacal pores (vs. C. 
auribalteatus 6, C. dumnuii 5 or 6, C. erythrops 9, C. huongsonensis 
6-8, C. kunyai 3, C puhuensis 5, and C. sonlaensis 8). In addition, 
C. phukhaensis sp. nov. is superficially similar to C. martini in 
having reticulations on the head but it differs from C. martini 
by being substantially larger in size (SVL: 77.8-82.8 mm vs. 
64.4-962 mm), having more precloacal pores (7 vs. 4), fewer 
supralabials (8 or 9 vs. 9-11). Furthermore, C. phukhaensis sp. 
nov. is distinguishable from C. zhenkangensis in having clear 
reticulations on the head and different body coloration in life 
(dark chocolate brown vs. light-brown); it also differs from 
C. zhenkangensis by being significantly smaller in size (SVL: 
77.8-82.8 mm vs. 87.4 mm), having fewer precloacal pores (7 
vs. 8), and possessing fewer transverse bands on the dorsum of 
the body (4 or 5 vs. 8 or 9). Cyrtodactylus phukhaensis sp. nov. 
differs from C. hekouensis in having reticulations on the head, 
dark body color, fewer precloacal pores (7 vs. 33-39), and fewer 
ventral scale rows (40-47 vs. 68-72) as well as by being smaller in 
size (SVL: 77.8-82.8 mm vs. 923 mm). 


3.4. Karyotype analysis The diploid chromosome numbers 
of C. phukhaensis sp. nov. were 2n = 40, and the fundamental 
number (NF) was 46. The karyotype of C. phukhaensis sp. nov. 
was composed of 2 metacentric, 4 submetacentric, and 34 
telocentric chromosomes (Figure 8). The karyotype formula was 
as follows: 2n (40) = 2m + 4sm + 34t. There was no sex difference 
between the male and female karyotypes. The ideogram of C. 
phukhaensis sp. nov. is shown in Figure 9 with the standard 
values given in Table 4. 


4. Discussion 
Cyrtodactylus phukhaensis sp. nov. was placed in the 


Cyrtodact ylus chauquangensis group because the type locality was 
in the area of distribution previously mentioned by Grismer 
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Figure 8 Conventional metaphase plates and karyotypes of male (A) and female (B) Cyrtodactylus phukhaensis sp. nov. (scale bar = 


10 pm). 


8 9 10 11 12 13 14 15 16 17 18 19 20 


Figure 9 Ideogram of Cyrtodactylus phukhaensis sp. nov. according to conventional staining. 


et al. (2021b), i.e., ranging from northern Thailand and Laos 
to northeastern Vietnam and Yunnan Province in China. 
However, the low genetic divergence between C. phukhaensis 
sp. nov. and C. wayakonei (5.2296) was surprising given that 
the distance between the type locality of these species is 
approximately 170 km, which is far greater than the distance 
between C. phukhensis sp. nov. and C. vilaphongi or Cyrtodactylus 
spelaeus (both approximately 125 km). This may be attributable 


to the adjacent mountain range between the Luang Prabang 
mountain range, on which Doi Phu Kha is located, and the 
mountains in Luang Nam Tha Province. This should be 
investigated further with the inclusion of information on the 
new Cyrtodact ylus spp. found on this mountain range. Thus, 
there is a need for field work on this endemic species in this 
region including the many hidden areas in northern Thailand. 
The chromosomal data on C. phukhaensis sp. nov. is shown 
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in Table 5 with data from previous reports on the genus 
from Thailand and Malaysia. The chromosome number for 
C. phukhaensis sp. nov., 2n = 40, was in agreement with those 
from other Cyrtodactylus species, which ranged from 32 to 
48 (King, 1987), showing that chromosomal diversification 
existed within the genus (Ota et al., 1992). In particular, the 
chromosomal numbers of C. phukhaensis sp. nov. were similar 
to those of C. kunyai and Cyrtodact ylus jarujini (2n = 40), one of 
the most frequent numbers found among Cyrtodact ylus spp., 
ie, 2n = 40 or 2n = 42 (Thongnetr et al, 2021). The obtained 
karyological characteristics of C. phukhaensis sp. nov. included 
many acrocentric chromosomes that gradually decreased in 
size (telocentric) and a few metacentric or submetacentric 


chromosomes. These features are consistent with the karyotype 
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evolution of other gekkonids and agree with the hypothesis of 
rearrangement from ancestral karyotypes via Robertsonian 
fissions, fusions, and pericentric inversions (Gorman, 1973). 
However, additional studies of cytogenetic approaches are 
required to clarify the chromosomal characteristics of not only 


C. phukhaensis sp. nov. but also other species from Cyrtodactylus. 
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Table 4 Mean length of short (p) and long (q) chromosome arm, total length of chromosome (total), centromeric index (CI) with standard 
deviation (SD), and relative length (RL) from metaphase of male and female Cyrtodactylus phukhaensis sp. nov. 


Chr. pair p+SD q+SD total + SD CI +SD RL Chr. type 
il 1.33 £0.14 2.01 + 0.20 3.34 € 0.29 39.74 + 2.64 0.1208 m 

2 0.88 + 0.13 1.63 + 0.24 2.51 € 0.33 35.11 x 3:65 0.0908 sm 

3 0.81 € 0.13 1521020 2.33 + 0.24 34.84 + 4.95 0.0842 sm 

4 0.00 + 0.00 2.25 + 0.23 2.25 € 023 0.00 + 0.00 0.0815 t 

5 0.00 + 0.00 2.08 + 0.25 2.08 + 0.25 0.00 + 0.00 0.0751 t 

6 0.00 + 0.00 1.92 + 0.24 1.92 + 0.24 0.00 + 0.00 0.0696 t 

z 0.00 + 0.00 1.68 + 0.26 1.68 € 0.26 0.00 + 0.00 0.0608 t 

8 0.00 + 0.00 1.57 € 0.23 1.57 € 0.23 0.00 + 0.00 0.0566 t 

9 0.00 + 0.00 1.42+0.18 1.42 £0.18 0.00 + 0.00 0.0515 t 

10 0.00 + 0.00 1.32 € 0.22 1:32£0.22 0.00 + 0.00 0.0476 t 

11 0.00 + 0.00 1.17 £ 0.18 1.17 £0.18 0.00 + 0.00 0.0423 t 

12 0.00 + 0.00 0.97 € 0.15 0.97 € 0.15 0.00 + 0.00 0.035 t 

13 0.00 + 0.00 0.87 + 0.14 0.87 + 0.14 0.00 + 0.00 0.0314 t 

14 0.00 + 0.00 0.79 £ 0.11 0.79 + 0.11 0.00 + 0.00 0.0284 t 

15 0.00 + 0.00 0.74+0.11 0.74 + 0.11 0.00 + 0.00 0.0269 t 

16 0.00 + 0.00 0.69 + 0.10 0.69 + 0.10 0.00 + 0.00 0.0249 t 

17 0.00 + 0.00 0.62 + 0.10 0.62 + 0.10 0.00 + 0.00 0.0223 t 

18 0.00 + 0.00 0.52 + 0.10 0.52 + 0.10 0.00 + 0.00 0.019 t 

19 0.00 + 0.00 0.47 + 0.08 0.47 + 0.08 0.00 + 0.00 0.0168 t 

20 0.00 + 0.00 0.40 + 0.08 0.40 + 0.08 0.00 + 0.00 0.0145 t 
Note: Chr. = chromosome; m = metacentric; sm = submetacentric; t = telocentric. 

Table 5 Review of karyotype studies for the genus Cyrtodactylus. 

Species 2n NF Karyotype formula NORs Location Reference 

C. consobrinus 48 50 2bi-arm+46t - Malaysia Ota et al. (1992) 

C. doisuthep 34 56 14m+6sm+2a+12t P9, 13 Thailand Thongnetr et al. (2021) 
C. interdigitalis 42 52 4m+2sm+4a+32t P12 Thailand Thongnetr et al. (2019a) 
C. jarujini 40 56 8m+4sm+4a+24t P13, 14 Thailand Thongnetr et al. (2021) 
C. kunyai 40 52 8m+4sm+6a+22t RIR Thailand Thongnetr et al. (2019a) 
C. phukhaensis sp. nov. 40 46 2m+4sm+34t - Thailand This study 

C. pubisulcus 42 44 2bi-arm--40t - Malaysia Ota et al. (1992) 

C. saiyok 42 42 42t P15 Thailand Thongnetr et al. (2019b) 


Note: 2n = diploid chromosome; NF = fundamental number; NOR = nucleolar organizer region; bi-armed = bi-armed chromosome; m = 


metacentric; sm = submetacentric; a = acrocentric; t = telocentric. 
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